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Abstract Achillides is a subgenus of the genus Papilio commonly called swallowtail butterflies, 
which inhabits East Asia including Japan. Achillides consists of many similar species, of which 
subspecies are distributed in many East Asian islands. Here, we analyzed phylogenetic relation- 
ships among these species and their subspecies with nucleotide sequences of the mitochodrial ND5 
gene. The analysis revealed that Papilio polyctor and P. bianor bianor are an identical species, and 
P. syfanius and P. maackii are also an identical species due to the identity of their nucleotide se- 
quences. These butterflies have been considered taxonomically different species based on their dis- 
tinct wing patterns. The analysis also revealed that Achillides are separated into two clades, which, 
however, do not agree with the grouping based on wing pattern similarity. One clade is composed 
of P. bianor, P. polyctor, P. hermeli, and P. dialis, while the other clade is composed of P. maackii, 
P. syfanius, P. krishna, P. hopponis, P. arcturus, P. karna and P. paris. Several pairs of species 
having a similar wing pattern are present in both clades. These findings may be evidence of paral- 
lel evolution, indicating that morphological convergence occurred between the species belonging to 
different clades after their ancestors had diverged. Speciation due to the wing pattern changes may 
have occurred in a relatively short period during the evolution of Achillides. Finally, we also dis- 
cuss the revision of classification of the P. bianor subspecies and P. fuscus based on DNA sequence 
analysis. 
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Introduction 


Molecular phylogeny of papilionid butterflies, commonly called swallowtails, has been ex- 
tensively studied recently (Sperling, 1993a, b; Caterino & Sperling, 1999; Yagi et al., 1999, 
Kato and Yagi, 2004; Zakharov et al., 2004), because their taxa and worldwide distributions 
were well elucidated. Achillides is a subgenus of the genus Papilio, which mainly inhabits 
east Asia. Achillides consists of at least twenty morphologically similar species all of 
which have black wings covered with beautiful bright green scales (Shimogori, 1997; Bauer 
& Frankenbach, 1998). Some species of Achillides are distributed in the far east Asian is- 
lands, where the species are classified into subspecies. 


One of the issues of the subgenus Achillides is the presence of taxonomically ambiguous 
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groups. The taxonomic relationships between Papilio polyctor and P. bianor and those be- 
tween P. syfanius and P. maackii have been disputed for a long time, because a cline of 
wing patterns between the species was found in some areas (Harada, 1992; Fujioka, 1997; 
Yoshimoto, 1998). Most taxonomists classified these groups to different species because 
their typical wing patterns are obviously distinct (Igarashi, 1979; Shimogori, 1997; Bauer & 
Frankenbach, 1998). 


The classification of P. fuscus is also ambiguous. The adult P. fuscus has the wing pattern 
that is a typical morphological feature of the subgenus Menelaides but its larva and pupa 
have morphological features of the subgenus Achillides. Igarashi (1979) classified P. fus- 
cus as belonging to Achillides. 


Another unresolved issue is the classification of P. bianor, which has a latitudinally wide 
distribution in east Asia. P. bianor is distributed from Russian Sakhalin Island, a subfrigid 
region, to Vietnam and Thailand, subtropical regions, through Korea, China, Japan and the 
Ryukyu Archipelago, temperate regions. P. bianor is classified into many subspecies (cf. 
Table 1), and some of them are suggested to be genetically distinct as species (Kawazoé 
and Wakabayashi, 1976). These observations prompted our interest in phylogenetic rela- 
tionships among the subspecies. 


Table 1. Taxa and locality of specimens of the subgenera Achillides and Menelaides used in this study”. 











species subspecies locality database accession number 
Achillides 

bianor bianor Fujian, China AB212916 
Sichuan, China AB212917 

Guangxi, China AB212918 

dehaanii Sakhalin, Russia AB212919 
Myohyangsan, North Korea AB212920 

Chunchon, South Korea AB212921 

Biei, Hokkaido, Japan AB212922 

Ishinomaki, Miyagi, Honshu, Japan AB212923 

Tohi-cho, Shizuoka, Honshu, Japan AB212924 

Shirahama-cho, Wakayama, Honshu, Japan AB212925 

Kyoto, Honshu, Japan ABO13156 

Kamiyama-cho, Tokushima, Shikoku, Japan AB212926 

Kasuya-cho, Fukuoka, Kyushu, Japan AB212927 

Tsushima Is., Japan AB212928 

tokaraensis Tokara-Nakanoshima Is., Japan AB212929 
Tokara-Akusekijima Is., Japan AB212931 

hachijonis Hachijo Is. 1, Japan AB212930 
Hachijo Is. 2, Japan AB223163 

Miyakejima Is. 1, Japan AB223164 

Miyakejima Is. 2, Japan AB223165 

Mikurajima Is., Japan AB223166 

thrasymedes — Nantou, Taiwan AB223167 
Taoyuan, Taiwan AB223168 

kotoensis Lanyu Is., Taiwan AB223169 
Huoshao Is., Taiwan AB223170 

okinawensis Iriomotejima Is., Japan AB231690 
Ishigakijima Is., Japan AB231691 

Taketomiima Is., Japan AB231692 

Yonagunijima Is. 1, Japan AB231693 

Yonagunijima Is. 2, Japan AB231694 
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Table 1 (continued). 


polyctor 


maackii 


syfanius 


hermeli 
dialis 
krishna 
hopponis 
arcturus 
paris 


karna 


fuscus 
Menelaides 
protenor 
helenus 
polytes 


amamiensis 


ryukyuensis 


polyctor 
stockleyi 
triumphator 


ganesa 
maackii 


tutanus 


albosyfanius 
kitawakii 


paris 
nakaharai 
hermosanus 
tamilana 
karna 
irauana 
carnatus 
palinurus 
peranthus 
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Amami-oshima Is. 1, Japan 
Amami-oshima Is. 2, Japan 
Tokunoshima Is. 1, Japan 
Tokunoshima Is. 2, Japan 
Okinawajima Is. 1, Japan 
Okinawajima Is. 2, Japan 
Kashmir, India 

West Thailand 

Laksao, Laos 

Nonhet, Laos 

Tibet, China 

Primorskyi, Russia 
Heilonjjiang, China 

Jilin, China 

Zhejian, China 

Sichuan, China 

Solaksan, South Korea 
Sakhalin, Russia 

Hokkaido, Japan 

Sendai, Miyagi, Honshu, Japan 
Yuzawa-cho, Niigata, Honshu, Japan 
Kyoto, Honshu, Japan 
Yamashiro-cho, Tokushima, Shikoku, Japan 
Kasuya-cho, Fukuoka, Kyushu, Japan 
Kikuchi, Kumamoto, Kyushu, Japan 
Tsushima Is., Japan 
Tanegashima Is., Japan 
Yakushima Is., Japan 

Yunnan, China 

East Tibet, China 

Mindoro Is., Philippine 
Nantou, Taiwan 

Nepal 

Nantou, Taiwan 

Sichuan, China 

Guangdong, China 

Nantou, Taiwan 

Nantou, Taiwan 

South India 

West Java, Indonesia 

Palawan, Philippine 

Sabah, Borneo, Malaysia 
Malaya, Malaysia 

Bali Is., Indonesia 

Celebes, Indonesia 


Kyoto, Honshu, Japan 
Kyoto, Honshu, Japan 
Iriomotejima Is., Japan 





AB239775 
AB239776 
AB239777 
AB239778 
AB013157 
AB239779 
AB239780 
AB239781 
AB239782 
AB239783 
AB239784 
AB239808 
AB239809 
AB239810 
AB239811 
AB239812 
AB239813 
AB239814 
AB239815 
AB239816 
AB239817 
AB239818 
AB239819 
AB013158 
AB239820 
AB239821 
AB239822 
AB239823 
AB239824 
AB239825 
AB239826 
AB239827 
AB239828 
AB239829 
AB239830 
AB239831 
AB239832 
AB239833 
AB239834 
AB239835 
AB239836 
AB239837 
AB239838 
AB239839 
AB239840 


AB013153 
AB013152 
AB013151 





» Revision of classification of P. bianor, P. polyctor, P. syfanius and P. fuscus is proposed in the text of 


this paper. 


d The names of P. bianor subspecies in this paper tentatively follow the nomenclature used in the books 
of Hidaka et al. (1984), Inomata and Matsumoto (1995) and Fujioka (1997). It should be noted that the 
other nomenclature used by Shir6zu (1966) and Kawazoé and Wakabayashi (1976) is also present. 
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To reveal taxonomic or phylogenetic relationships among these rather complicated species 
and their subspecies of the subgenus Achillides, nucleotide sequences of mitochondrial 
DNA should be very useful. The mitochondrial NADH dehydrogenase subunit 5 (ND5) 
gene is one of the fastest-evolving genes (Garesse, 1988) and has been used for the molecu- 
lar phylogenetic analysis of insects (Su, Tominaga et al., 1996; Su, Ohama et al., 1996; Su 
et al., 1998). We previously showed that NDS sequences could be used to reveal the phylo- 
genetic relationships among lower categories of Lepidoptera (Yagi et al., 1999). 


In this paper, we show phylogenetic trees of most species and their subspecies of Achillides 
based on NDS nucleotide sequences, and discuss the relationship between phylogenetic di- 
vergence and wing pattern evolution. We also show NDS phylogenetic trees of P. bianor 
collected at various localities in east Asia, and discuss the classification and biogeography 
of the species relative to the paleogeography of the Ryukyus and Japan, the far east Asian 
islands. 


Materials and Methods 


Most species of the subgenus Achillides as 1-5-year-old dried specimens were kindly sup- 
plied by various Japanese entomologists (Table 1). Most butterflies collected by the au- 
thors, mainly in Japan, had been stored in ethanol until use. DNA was extracted from mus- 
cles of the thorax and legs by the phenol-chloroform extraction method described previous- 
ly (Yagi et al., 1999). In case of the mounted specimens, a few legs were removed for DNA 
extraction. 


A section of the mitochondrial NDS gene (810 bases) was amplified from the DNA by poly- 
merase chain reaction (PCR) using primers V1 and A1, also described previously (Yagi et 
al., 1999). Nucleotide sequences of the PCR products were directly determined with the 
BigDyeTerminator Cycle Sequencing FS Kit by an automated DNA sequencer 373A or 377 
(Applied Biosystems). The primers V1 and Al were usually used for the nucleotide se- 
quencing, and internal primers A3 (5’-TTCGAATTTAGCTTTATGTGG-3’) and C2 (5’-AT- 
CYTTWGAATAAAAYCCAGC-3’) were also used when V1 and A1 did not produce clear 
results. 


NDS nucleotide sequences (789 bp) of the samples were edited and aligned using the 
YooEdit software (version 1.63). Phylogenetic trees were constructed with the neighbor- 
joining (NJ) method and the unweighted pair-group method using an arithmetic average 
(UPGMA) by the PHYLIP software package version 3.572c (Felsenstein, 1993), or with 
the maximum parsimony (MP) method by the PAUP software version 4.0 beta (Swofford, 
1998). For NJ and MP methods, P. memnon, P. protenor and P. helenus belonging to an- 
other subgenus Menelaides were taken as outgroup species. The confidence level of 
branching in the phylogenetic trees (NJ and UPGMA) was evaluated with the bootstrap test 
based on 100 resamplings (Felsenstein, 1985). The nucleotide sequences of the ND5 gene 
of these butterflies were registered in the DDBJ/GenBank/EMBL database. 


Results 


Alignment of the NDS nucleotide sequences from Achillides species and their subspecies 
revealed that DNA length variation was not present and the maximum pairwise base differ- 
ences were 6.2%, suggesting that NDS sequences are valid for the phylogenetic analysis of 
Achillides. The phylogenetic analyses of Achillides species and their subspecies were car- 
ried out with the NJ, MP and UPGMA methods. As all these methods gave comparative 
phylogenetic trees, only the NJ tree is shown in Fig. 1. Fig. 1 indicates that P. peranthus 





NII-Electronic Library Service 


The Lepidopterological Society of Japan 


Molecular Phylogenetic Analysis of Subgenus Achillides 141 


51 


100 


77 


0.01 


79 
@ 73 =f. bianor kotoensis Lan-Yu Is., Taiwan 


P. syfanius albasyfanius Yunnan, China 
P. maackii tutanus Sakhalin, Russia 
P. maackii tutanus Yamashiro-cho, Japan 
P. maackii tutanus Yakushima |s., Japan 
P. maackii tutanus Tanegashima Is., Japan 
P. maackii tutanus Tsushima is., Japan 
P. maackii maackii Jilin, China 
P. maackii maackii Primorskiy, Russia 
P. maackii maackii Zhejian, China 
P. maackii tutanus Yuzawa-cho, Japan 
P. maackii maackii Solaksan, S. Korea 
P. maackii maackii Sichuan, China 
P. maackii maackii Heilongjiang, China 
P. maackii tutanus Fukuoka, Japan 
614 PF. maackii tutanus Kikuchi, Japan 
100 P. maackii tutanus Sendai, Japan 
P. maackii tutanus Kyoto, Japan 


65 P. maackii tutanus Hokkaido, Japan 
74 P. syfanius kitawakii E. Tibet, China 
P. Krishna Nepal 
92 P. hopponis Nantou, Taiwan 
P. arcturus Sichuan, China 
66 5a 100 „2. karna carnatus Borneo, Malaysia 


P, karna irauana Palawan, Philippine 


D P. karna karna W. Java, indonesia 


P. paris hermosanus Nantou, Taiwan 
100 Í 99 p £. paris paris Guangdong, China 
P. paris nakaharai Nantou, Taiwan 
P, paris tamifana S. india 
P. bianor dehaani Sakhalin, Russia 
P. bianar dehaani Kyoto, Japan 
100 P. bianor dehaanii Chunchon, S. Korea 
P. bianor tokaraensis Akuseki-jima Is., Japan 
P. bianor hachijonis Mikura-jima is., Japan 
P. bianor amamiensis Amami-ashima Is., Japan 
P. bianor amamiensis Tokunoshima is., Japan 
P. bianor ryukyuensis Okinawa-jima Is., Japan 
P. hermeli Mindoro ts., Philippine 
P. dialis Nantou, Taiwan 
P. bianor bianor Fujian, China 


100 


P polyctor triumphator Laos 
P polyctor stockleyiW. Thailand 
P. bianor thrasymedes Nantou, Taiwan 
100 =”. bianor okinawensis \riomote-jima Is., Japan 
£. bianor okinawensis Yonaguni-jima Is., Japan 
P., palinurus Malaya, Malaysia 
P. peranthus Bali, Indonesia 
P. polytes 
P. protenor 
P. fuscus Celebes, indonesia 
P. helenus 


Fig. 1. The phylogenetic tree of Achillides butterflies constructed with the NJ method using the 
PHYLIP software package. P. helenus belonging to another subgenus Menelaides was 
taken as an outgroup. The confidence level of branching in the phylogenetic trees was eval- 
uated with the bootstrap test based on 100 resamplings, and the values higher than 40 are 
shown at branching points of the tree. In Achillides, P. peranthus and P. palinulus diverge 
early, and the so-called real Achillides group then diverges to two groups, which are named 


clades 1 and 2 in the tree. 


and P. palinulus diverged earlier than the so-called genuine Achillides group in the sub- 
genus Achillides. The genuine Achillides then diverged to two groups (clades 1 and 2 in 
Fig. 1). Clade 1 consists of P. dialis, P. hermeli, P. bianor and P. polyctor, while clade 2 
consists of P. paris, P. karna, P. arcturus, P. hopponis, P. krishna, P. maackii and P. syfa- 
nius. P. fuscus was not included in Achillides but in the Menelaides group. 


To clarify the phylogenetic relationships among the subspecies of P. bianor, butterflies col- 
lected at various localities were analyzed with the NJ and UPGMA methods. These meth- 
ods gave phylogenetic trees with a similar topology except for the divergences with low 
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a Sakhalin, Russia 
My ohyangsan, N. Korea 
Chunchon, S. Korea 
Biei, Hokkaido, Japan 
Ishinomaki, Honshu, Japan 
Tohi, Honshu, Japan P. bianor dehaanii 
Shirahama, Honshu, Japan 
Kyoto, Honshu, Japan 
Kamiyama, Shikoku, Japan 
Kasuya, Kyushu, Japan 
Tsushima Is. Japan 
Tokara-nakanoshima Is., Japan 
Tokara-akusekijima Is., Japan 
Hachijo-jima Is. 1, Japan 
100 Hachijo-jima Is. 2, Japan 
Miyake-jima is. 1, Japan | ZP. bianor hachijonis 
Miyake-jima is. 2, Japan 
Mikura-jima Is., Japan 
47 Tokunoshima Is. 1, Japan 
70 Tokunoshima Is. 2, Japan 
100 Amami-oshima is. 1, Japan 
Amami-oshima Is. 2, Japan 
69 Okinawa-jima Is. 1, Japan 
Okinawa-jima Is. 2, Japan 
Fujian, China 
Sichuan, China 
Laksao, Laos ~ P. polyctor triumphator 
42 W. Thailand — P. potyctor stockleyi 
Kashmir, India === P, polyctor polyctor 
61 FHuoshao Is., Taiwan 
Lanyu Is., Taiwan 
60 Nonhet, Laos —— PF, polyctor triumphator 
Nantou, Taiwan 
100 Taoyuan, Taiwan 
Guangxi, China — P. bianor bianor 
Tibet, China —— P, polyctor ganesa 
52 lriomote-jima Is., Japan 
Ishigaki-jima Is., Japan 
pos Taketomi-jima iS., Japan | Z. bianor okinawensis 
Yonaguni-jima Is. 2, Japan 
Yonaqguni-jima Is. 1, Japan 
P. dialis 
80 P. maackii 
99 P. krishna 
P. hopponis 


| P. bianor tokaraensis 


P. bianor amamiensis 


| P. bianor ryukyuensis 


| P. bianor bianor 


| 2. bianor kotoensis 


| P. bianor thrasymedes 


0.01 


Fig. 2. The phylogenetic trees of P. bianor and P. polyctor collected at various regions. The trees 


were constructed with the NJ (a) and UPGMA (b) methods by the PHYLIP software pack- 
age. P. hopponis, P. krishna and P. maackii were used as outgroup species in the NJ tree. 
The confidence level of branching in the phylogenetic trees was evaluated with the boot- 
strap test based on 100 resamplings, and the values higher than 40 are shown at branching 
points of the trees. 


bootstrap values (Fig. 2). As the evolution rate of individuals would be almost constant 
within a species, the UPGMA tree is used further to discuss the biogeography of P. bianor. 


Discussion 


Revision of classification by DNA sequences 


Figs 1 and 2 unexpectedly show that nucleotide sequences of P. polyctor including three 
subspecies polyctor, stockleyi and triumphator coincide with those of P. bianor bianor, and 
the sequences of P. syfanius coincide with those of P. maackii. P. polyctor has a large sky- 
blue spot on the apical area of the hindwings but P. bianor does not (Fig. 3). P. syfanius 
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Sakhalin, Russia 
b My ohyangsan, N. Korea 
Chunchon, S. Korea 
Biei, Hokkaido, Japan 
ishinomaki, Honshu, Japan 
Tohi, Honshu, Japan P. bianor 
Shirahama, Honshu, Japan dehaanii 
Kyoto, Honshu, Japan 
Kamiyama, Shikoku, Japan 
Kasuya, Kyushu, Japan 
Tsushima Is. Japan 
Tokara-nakanoshima Is., Japan} 7. bianor 
Tokara-akusekijima Is., Japan tokaraensis 
Hachijo-jima Is. 1, Japan 
Hachijo-jima Is. 2, Japan 


= Miyake-jima Is. 1, Japan | 9 S/anor hachyonis 
Miyake-jima Is. 2, Japan 
Mikura-jima is., Japan 
Fujian, China | P. bianor bianor 
Sichuan, China 
Laksao, Laos — P. palyctor triumphator 
W. Thailand ~m P, polyctor stockleyi 

49 Kashmir, india ~~ F. polyctor polyctor 

Huoshac:Is., Taiwan | 2. bianor kotoensis 
Lanyu Is., Taiwan 

44 Guangxi, China — P. bianor bianor 
Nonhet, Laos wm P. polyctor triumphator 

= Tibet, China — P. polyctor ganesa 


Nantou, Taiwan 


Taoyuan, Taiwan | P. bianor thrasymedes 


93 55 Iriomote-jima Is., Japan 
Ishigaki-jima Is., Japan . 
57 
100 g | Taketomi-jima Is., Japan P. bianor 
okinawensis 


Yonaguni-jima Is. 2, Japan 
Yonaguni-jima Is. 1, Japan 
Amami-oshima Is. 1, Japan 

7S fI Amami-oshima Is. 2, Japan | 2, bianor 

100 Tokunoshima Is. 1, Japan amamiensis 

Tokunoshima Is. 2, Japan 
Okinawa-jima Is. 1, Japan P, bianor 

70 bm Okinawa-jima Is. 2, Japan ryukyuensis 
P. dialis 
P. hopponis 
P. maackii 


P. krishna 
0.01 


Fig. 2 (continued). 


has a large white spot on the postmedian area of the hindwings but P. maackii does not 
(Fig. 3). Therefore, these butterflies have been considered taxonomically different species 
due to their distinct wing patterns. Here, our molecular phylogenetic analysis revealed that 
P. polyctor and P. bianor bianor should be an identical species, and P. syfanius and P. 
maackii should also be an identical species. These findings are reinforced by the evidence 
that distribution of the two species is roughly allopatric in the Eurasian Continent and a 
cline of their wing patterns has been observed in the small area where their distributions 
overlap (Fujioka, 1997; Shimogori, 1997). 


P. fuscus was found to be a species of the subgenus Menelaides from mitochondrial DNA 
sequences. P. fuscus has been included in the subgenus Achillides since the description by 
Igarashi (1979), because the morphological features of the larva and pupa are similar to 
those of Achillides species. Our result shows that the wing pattern of P. fuscus reflects its 
phylogenetic features. 
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d e 


Fig. 3. Wing pattern convergence between phylogenetically distinct species. a: P. bianor bianor 
(Sichuan, China); b: P. polyctor polyctor (North India); c: P. syfanius albosyfanius 
(Yunnan, China); d: P. maackii maackii (Sichan, China); e: P. paris paris (North India); f: 
Parides nevilli (Yunnan, China). Upper and lower individuals in each vertical row inhabit 
the same region and have similar wing patterns but are phylogenetically distinct. P. bianor 
bianor (a) and P. polyctor polyctor (b), and also P. syfanius (c) and P. maackii (d) should be 
identical species, although they have quite different wing patterns. 


In Fig. 1, bootstrap values higher than 40 were shown at branching points of the phyloge- 
netic tree. It should be noted that the branching with the value lower than approximately 60 
is unreliable. For instance, connection of P. paris hermosanus to the P. karna group, and 
inclusion of P. hermeli and P. dialis in the P. bianor group are not trustworthy. Other di- 
vergences at the points with low bootstrap values in the phylogenetic trees also could not be 
discussed. 


P. bianor was divided into four groups (Figs 1 and 2), whose pair-wise Kimura’s two-pa- 
rameter distances were larger than 0.019. The first group includes ssp. dehaanii, tokaraen- 
sis, hachijonis, which inhabits the Japan Archipelago and the Korean Peninsula. The sec- 
ond group includes ssp. amamiensis and ryukyuensis, which inhabit the middle Ryukyu 
Islands. The third group includes ssp. bianor, thrasymedes, kotoensis and P. polyctor, 
which inhabit continental China and the Taiwan Islands. The fourth group consists of ssp. 
okinawensis, which inhabit the Yaeyama Islands in south Ryukyu. These four groups 
should not be regarded as subspecies but four independent species, because their pair-wise 
genetic distances are as large as the distances between other distinct species, for instance, 
between P. paris and P. karna, and between P. hopponis and P. arcturus (Fig 1). In addi- 
tion, the divergence of four P. bianor groups took place at almost the same time as those of 
other species P. dialis and P. hermeli (Fig. 1). This consideration is supported by failure of 
normal development of pair-wise hybrids of these four group individuals (Ae, 1995; 
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Kawazoé, pers. comm.). The divergence order of the four groups was not clarified in this 
study, because the bootstrap values were low, and therefore, the order is variable among the 
analytical methods (Fig. 2). 


Parallel evolution and convergence of wing pattern 


Similarity of butterfly wing patterns contradicts the mitochondrial DNA phylogeny in the 
subgenus Achillides. Pairs of species with a similar wing pattern are present between 
clades 1 and 2, for example; P. bianor bianor and P. maackii maackii, P. polyctor polyctor 
and P. paris paris (Fig. 1). Fig. 3 (a and d; b and e) shows pairs of species of clades 1 and 
2 with a similar wing pattern. These pairs of species are distributed in identical or neigh- 
boring areas (Fujioka, 1997; Shimogori, 1997). The phylogenetic tree (Fig. 1) indicates that 
phylogenetic divergences occurred earlier than wing pattern changes, and the wing patterns 
of the different species then converged independent of the phylogeny. The species inhabit- 
ing a geographically isolated region may be liable to develop a different wing pattern from 
the original one, and some different species inhabiting the same region or a similar environ- 
ment may evolve a similar pattern on their wings. These findings are similar to the phe- 
nomenon, parallel evolution, found in races of South American long wing butterflies 
Heliconius erato (Brower, 1994a, b), Japanese ground beetles Carabus (Su, Tominaga et 
al., 1996; Su, Ohama et al., 1996) and African lake fish cichlid (Meyer et al., 1990), sug- 
gesting that morphologically similar species inhabiting different areas are phylogenetically 
distinct, and that morphologically different species inhabiting the same area are phylogenet- 
ically proximate. In any case, these studies revealed that morphological changes of species 
occurred independently of their phylogenetic divergences. 


NDS evolutionary rate and biogeography of P. bianor 


As the UPGMA tree is thought to be correlated with the past phylogenetic divergence time 
in a lower taxonomic group, it would be reasonable to hypothesize that the branching points 
of the four P. bianor groups (Fig. 2) were settled at the time when the groups were isolated 
in the Japan and Ryukyu Archipelagoes. According to the paleogeography of the Ryukyu 
Archipelago, the Ryukyu and the Sakishima Super-islands formed by the end of the 
Pliocene (1.65 MYA), as a result of subsidence of the east periphery of the Eurasian 
Continent during the Pliocene (1.65-5.3 MYA) (Kizaki & Oshiro, 1977; Ujiié, 1990; Ota, 
1998). The Ryukyu and Sakishima Super-islands separated far later to the present Okinawa 
and Amami Islands, and the present Yaeyama and Miyako Islands, respectively (Kizaki & 
Oshiro, 1977; Ota, 1998). From this assumption, Kimura’s two-parameter distance 0.01D 
would be estimated to be about 0.7-0.8 million years (MY). This rough estimation can also 
reasonably explain that the divergence points of the continental and Taiwan populations, of 
the Iriomote and Yonaguni populations, and of the Okinawajima and Amami-Oshima popu- 
lations in the UPGMA tree are about 0.17-0.35 MYA, which roughly corresponds to the 
time the present Ryukyu Islands formed in the late Pleistocene (Kizaki & Oshiro, 1977; 
Ujiié, 1990; Ota, 1998). This NDS evolutionary rate of P. bianor agrees with that of 
Parnassius butterflies (0.01D=0.75 MY) (Yagi et al., 2001). 


The distinct wing pattern difference between P. polyctor and P. bianor bianor, and between 
P. syfanius and P. maackii without NDS sequence change indicates that the wing pattern 
change occurred in the period while one base change mutation had been fixed in the NDS 
gene (shorter than 0.07 MY). The wing pattern of P. polyctor polyctor resembles that of 
the phylogenetically distinct P. paris paris, and P. syfanius may mimic the poisonous papil- 
jonid species Parides nevilli, Par. hedistus and Par. polyeuctes (tribe Troidini) which occur 
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sympatrically (Figs 3c and f). These results indicate that wing pattern convergence in 
Achillides may occur in a relatively short period. 
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H Ea 


S hay FY PT DNAMPIC ko THOPIA TOR (TTY a ET 
F a YE) ORK EWE, EKORRAR ORZ + eX Nill ERR) 


ARERR UT ICEIDAT ANTI AT TOE AIO ARREST bY PUY DNADNDS 
ESTÉT ko CHL. ZORE, WI ATTO\GRBILZOO FA V—TIKE lDN 
CEPR MSIE OD Va TINH AT FIG ORS FIT AS PEDGAS EIN IAD? 
BIAT FRBSSEUNA: RIOT MATIC S VRP I AT TIIN SFR ATION GRADY 
VDT PIN, AAD LRIT TIN RIRTITN, MY EYL TAN, ANVINMYU EVI TFARSENS. 
BT WV — SP ORE ICSE OSE DSH SZ DU Cloak <<, PV PHIL § OPORMAW BOA 
ABDUPAETA. PERISHED ODAIATANESV VATA ANALY FOZ SYY 
TRNEMY EVI TAR ECHSZ. CDE LIL2ZFWV— FORINT V—- SHORE A LU CHO 
BEM OUD COKE ZA ON, FATHOM Arh LUR. RAEO IAT}, 
DDTPM-TIZKESKBONACERBD ROK. TRedbSB ls bh A FIBDVILO BATE - Y\)) 
Y- RAPE, BIBRA: W25 MBE, BISBEE, PU HEKER - 
GAB BE CHS. CN6547V—TFODNABINOEB, FAVHCHA LL TOHYRIEEKS 
VIER, DIAS PNZRE LZ YYZ TTOODNARSI, TVVAFATRANEVIAIATST 
ID DNA BCS CMI EA EHR, CNSILEN EN -FECHS OL HR IETS. KR 
YVIAEVAET AP NNDB EOE R SOADATTO GB LANAT EEDA, VUFE 
TT IHIBETNACHALCLEMPDPOR. 
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